Equilibrium long-ranged charge correlations at the interface 
between media coupled to the electromagnetic radiation 
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We continue studying long-ranged quantum correlations of surface charge densities on the interface 
between two media of distinct dielectric functions which are in thermal equilibrium with the radiated 
electromagnetic field. Two regimes are considered: the non-retarded one with the speed of light 
c taken to be infinitely large and the retarded one with finite value of c. The analysis is based 
on results obtained by using fluctuational electrodynamics in [L. Samaj and B. Jancovici, Phys. 
Rev. E 78, 051119 (2008)]. Using an integration method in the complex plane and the general 
analytic properties of dielectric functions in the frequency upper half-plane, we derive explicit forms 
of prefactors to the long-range decay of the surface charge correlation functions for all possible 
media (conductor, dielectric, vacuum) configurations. The main result is that the time-dependent 
quantum prefactor in the retarded regime takes its static classical form, for any temperature. 

PACS numbers; 05.30.-d, 52.40.Db, 73.20.Mf, 05.40.-a 



I. INTRODUCTION 

In this paper, we continue studying long-ranged quan- 
tum correlations of surface charge densities on the inter- 
face between two distinct media, initiated in Refs. 
The model, formulated in the three-dimensional (3D) 
Cartesian space of coordinates (x, y, z), is inhomogeneous 
say along the first coordinate x (see Fig. 1). The two 
semi-infinite media with the frequency-dependent dielec- 
tric functions ei{uj) and e2(w) (the magnetic permeability 

— I) are localized in the half-spaces a; > and a; < 0, 
respectively. The interface is the plane x = 0, a point 
on the interface is R = (0, y, z). The different electric 
properties of the media give rise to a surface charge den- 
sity which must be understood as being the microscopic 
volume charge density integrated on some microscopic 
depth. It is related to the discontinuity of the x com- 
ponent of the electric field on the interface. Denoting by 
cr(t, R) the surface charge density at time t and at a point 
R, the (symmetrized) two-point correlation function, at 
times different by t, reads 



eico) 



e,(co) 



5(t,R) = i(a(t,R)a(0,O) 



^a(0,O)a(t,R))T, (1) 







FIG. 1: Two semi-infinite media characterized by dielectric 
functions ei(a;) and e2{uj). 



two cases are treated separately. 

The two media configuration studied so far was re- 
stricted to a conductor, localized say in the half-space 
X > with the dielectric function eiito) = e{Lo), in con- 
tact with vacuum of the dielectric constant e2{'-^) = 1 (the 
vacuum part of the space is equivalent to a hard wall im- 
penetrable to charged particles forming the conductor). 
In some theoretical studies, the dielectric function is a_ 
proximated by a simple one-resonance Drude formula 



X 



where (• • • represents a truncated statistical average at 
the inverse temperature /? (we exclude the case of zero 
temperature, P —^ oo). We are interested in the behavior 
of the correlation function ([1]) at distances on the inter- 
face i? = |R| large compared to the microscopic length 
scales (like the particle correlation function). The static 
case of zero time difference t = between the two dis- 
tinct points is simpler than the one with t and so the 



e(cj) = 1 + 



* Unite Mixte de Recherche No 8627-CNRS 

'Electronic address: Bernard. Jancovici@th.u-psud.frl 

tOn leave from Institute of Physics, Slovak Academy of Sciences, 

Bratislava; Electronic address: ,Ladislav.Samaj@savba.sk| 



lu'q — uj{lu + irj) ' 



(2) 



where ujp is the plasma frequency, rj the dissipation con- 
stant and luq the oscillation frequency of harmonically 
bound charges: loq = for conductors and luq for 
dielectrics. In [l| and 0, we apphed ^ to the jellium 
model (sometimes called the one-component plasma), i.e. 
a system of pointlike particles of charge e, mass m, and 
bulk number density n, immersed in a uniform neutral- 
izing background of charge density —en. The dynamical 
properties of the jellium have a special feature: There is 
no viscous damping of the long- wavelength plasma oscil- 
lations for identically charged particles, so that 0+ 
in (21). The frequencies of non-retarded nondispersive 



2 



long-wavelength collective modes, namely Up of the bulk 
plasmons and ajg of the surface plasmons, are given by 



\ m 



1/2 



71' 



(3) 



For other conductors and dielectrics, we shall use only 
general properties of e(a;), without restricting ourselves 
to the Drude model 

The problem of a conductor in contact with vacuum 
was studied in the past, with an increasing level of phys- 
ical complexity and their range of validity. The prob- 
lem can be treated as classical or quantum, in the non- 
retarded or retarded regime. In all cases, the asymptotic 
large-distance behavior of the surface charge correlation 
function ([1]) exhibits a long-ranged tail of type 



R 



(4) 



where the form of the prefactor function h[t) depends on 
the physical theory which is used. It is useful to introduce 
the Fourier transform 



S'(t,q)= j (fRexp{iq-Il)S{t,R), 



(5) 



with q — {qy,qz) being a 2D wave vector. Since, in the 
sense of distributions, the 2D Fourier transform of 1/R^ 
is —2iTq, a result equivalent to Q is that pS{t,q) has a 
kink singularity at q = 0, behaving like 



l3S{t,q) 2TTh{t)q, 



0. 



(6) 



• Classical non-retarded regime: Let the conduc- 
tor be modelled by a classical Coulomb fluid composed 
of charged particles with the instantaneous Coulomb in- 
teractions. By a microscopic analysis [4], the long-range 
decay of the static surface correlation function was found 
such that 



1 



(7) 



where the subscript "cl" means "classical" and the upper- 
script "nr" means "non-retarded" (i.e. considered with- 
out relativistic effects associated with the finiteness of 
the speed of light c) . The same result has been obtained 
later Q by simple macroscopic arguments based on a 
combination of the linear response theory and the elec- 
trostatic method of images. Note the universal form of 
/ici(O), independent of the composition of the Coulomb 
fluid. 

• Quantum non-retarded regime: The extension 
of the static result ([7]) to a quantum Coulomb fluid, mod- 
elled by the jellium, was accomplished in Ref. [6]. The 
absence of damping was crucial in the treatment using 
long-wavelength collective modes, the bulk and surface 
plasmons with frequencies given in ([3|). The Maxwell 
equations, obeyed by the plasmons, were considered in 
the non-retarded (non-relativistic) regime with the speed 



of light c taken to be infinitely large, c = oo, ignoring in 
this way magnetic forces acting on the charged particles. 
The obtained time-dependent result has the nonuniversal 
form [1,0, i 

h^^^it) = — 



X [2g{tJs) cos{uJst) — g{ujp) cos{ujpt)] , (8a) 

(8b) 



= — coth ( — 



where the subscript "qu" means "quantum" . Accord- 
ing to the correspondence principle, a quantum system 
admits the classical statistical description in the high- 
temperature limit f3h —^ 0. In this limit, the function 
g{LLj) = 1 for any to and the quantum formula (|8ap re- 
duces to the classical non-retarded one 



,(nr) 
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'8^ 



[2 cos(ti>st) — cos{ujpt)] . 



(9) 



For t = 0, we recover the classical static formula ([T]). 

• Quantum retarded regime: In the previous paper 
[l| , we studied the surface charge correlations taking into 
account retardation (c is assumed finite) and the quan- 
tum nature of both the jellium and the radiated electro- 
magnetic (EM) field which are in thermal equilibrium. In 
other words, the quantum particles are fully coupled to 
both electric and magnetic parts of the radiated EM field. 
By using Rytov's fluctuational electrodynamics [9|, Hfl] we 
showed that there are two regions of distances R on the 
interface: the intermediate one given by the inequalities 
Aph ^ i? <C c/ujp (Aph oc Phc stands for the thermal de 
Broglie wavelength of photon), where the non-retarded 
result ([8a| . ([Fb]) applies, and the strictly asymptotic one 
given by the inequality c/up <C i?, where a retarded re- 
sult applies. After long calculations in [l|, a bit shortened 
through the alternative method of , a very simple form 
of the retarded result was found 



1 

8^' 



(10) 



Here, the upperscript "r" means "retarded". We see that, 
for any temperature /? and time t, the inclusion of retar- 
dation effects causes the prefactor function h{t) to take 
its universal static classical form ([7|), independent of h 
and c. The formula ([TU|) does not change in the classical 
limit f3h — > 0, i.e. 



/iW(i) ^ 1,^. 

cl V ) 
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(11) 



The presence or the absence of magnetic fields is impor- 
tant in two-point classical statistical averages, taken at 
two different times [T]|. Therefore, it is not surprising 
that the classical retarded formula pTjl and the classical 
non-retarded one © do not coincide with one another 

for nonzero time differences, ^^i (^) ''■ci (^) ^ 7^ 0- 
On the other hand, for t = we have 



(0). 



(12) 
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This equality is in agreement with the Bohr- van Leeuwen 
theorem [ij, [3] about an effective ehmination of mag- 
netic degrees of freedom from statistical averages (with 
zero time differences among the fixed points in the coordi- 
nate part of the configuration space) of classical systems; 
for a detailed treatment of this subject, see Ref. 

In the previous papers P, we used the fact, spe- 
cial to the jellium model, that there is no damping for 
small wave numbers. The question whether the crucial 
formula (jTU]) is still valid for a conductor with dissipa- 
tion was left as an open problem. There are still many 
other unsolved physical situations which deserve atten- 
tion. What happens in the case of a general dielectric 
in contact with vacuum? Filling the vacuum region by 
a material medium, other types of contacts are possible, 
like conductor - conductor, dielectric - dielectric, conduc- 
tor - dielectric. Another kind of problem is the algebraic 
complicacy connected with the derivation of the result 
(llOP for the jellium with the relatively simple form of the 
dielectric constant. Does there exists a simple method for 
evaluating the large-distance asymptotics of the surface 
charge density correlation function which is applicable to 
the jellium as well as to other more complicated systems? 
All asked questions are answered in the present paper. 

The generalization of the formalism to contacts be- 
tween all kinds of materials, defined by their dielectric 
functions, has already been done in Ref. The true 
problem is the mathematical handling of the final (non- 
retarded or retarded) formula for the Fourier transform 
f3S(t,q), written as an integral over real frequencies, to 
deduce the small-g behavior Here, we accomplish the 
task first by extending the integration over real frequen- 
cies to a contour integration in the complex frequency 
plane and then using integration techniques in the com- 
plex plane together with the known analytic properties of 
dielectric functions in the frequency upper half-plane. In 
short, the time-dependent retarded results maintain the 
simplicity of the jellium formula (jlOp and involve only di- 
electric functions of media in contact at zero frequency. 
The non-retarded results are complicated and available 
only, in general, as infinite series over Matsubara frequen- 
cies. 

The article is organized as follows. Sec. II is a gener- 
alization of the classical static result ([7]) for conductor to 
a dielectric. The result will serve us as a check of more 
general calculations. Sec. Ill summarizes briefly general 
analytic properties of dielectric functions in the complex 
frequency upper half-plane which are necessary for the 
derivation of our basic results. Sec. IV applies the new 
method for the calculation of /3S'(0, q). Sec. V does the 
same for pS{t,q). Sec. VI is a Conclusion. 



II. STATIC CORRELATIONS FOR CLASSICAL 
DIELECTRICS 

It has been known for a long time Q that the classical 
surface charge correlations on a conductor, made of par- 



ticles interacting through the Coulomb law and bounded 
by a plane wall, are long-ranged. We prefer to use the 
macroscopic language |5] . If R and R' are two points on 
the wall, the classical correlation, for distances |R— R'| 
large compared to the microscopic scale, behaves like 



/3S(R - R') = /3(cr(R)cr(R')) - - 



1 



87r2|R-R' 



'13 



(13) 



(we assume that the conductor is uncharged, (ct(R)) = 
0). 

A generalization of for the case of a dielectric of 
static dielectric constant e(0) — eg bounded by a plane 
can be found by the same method which has been used in 
cr(R) is related to the discontinuity of the x compo- 
nent of the electric field on the wall. We call i?™*''"''(R) 
the limit of that field component as R is from the out- 
side (inside) of the dielectric. The correlation of the sur- 
face charge densities is 



(a(R)a(R')) 



1 



f&"(R) 



E° 



(R) 



(47r)2 

x[4"(R')-^^r(R')]>- (14) 



If a test infinitesimal charge is introduced at r = 
(x,0, 0), a; > in the dielectric, the electric potential 
created by the dielectric at a point r' = (x' , y' , z'), x' > 
in the dielectric is given by the method of images as 



1 



9(1 -eo) 



|r'-r| (i + eo)|r'-r* 



|r' - r| 
(15) 

where r* = {—x, 0, 0) is the image of r. The last term in 
is the potential created by q, which should not be 
included in the potential created by the dielectric. The 
linear response theory relates the response (ITSt to the 
unperturbed correlation function between the additional 
Hamiltonian q4>{v) and (/>(r'), giving 



W)), = -/3g(</)(r)<^(r')). 



(16) 



Since the x component of the electric field is Ex{r) — 
— {d/dx)(f){r), using ([15]) and p6|) we can obtain the cor- 
relation of the X component of the electrical fields inside 
the dielectric at the wall 



/3(£;i"(R)i?i"(R')) = 



1 



£0 



2 



1 



1 



|R-R'P- 



(17) 

Similar calculations give the correlation outside the di- 
electric 



and the cross-correlation 

(R)£;°"'(R')) - - 



l-eo 



1 



1 + eo |R-R'|3 



1 



1 



1 + eo |R-R'|3' 



(18) 



(19) 
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Using ini), (UHl) and ^ in ((14]), we obtain 



/3(a(R)a(R')) 



1 / 1 



1 



1 



87r2 Veo 1 + eo/ |R-R'|3' 

(20) 

which is the wanted generahzation for a dielectric. With 
regard to the definition we have the nonuniversal 
result 



/lel(O) 



1/1 4 
1 



87r2 V^o 1 + eo 



(21) 



The value of e(0) is infinite for any kind of conductor and 
we retrieve (US]) from or Q from jH]). This fact 

explains the universality of the static hd for conductors. 



III. ANALYTIC PROPERTIES OF DIELECTRIC 
FUNCTIONS 

Analytic properties of dielectric functions are described 
in many textbooks [3, H, [iBl • They apply to an arbitrary 
dielectric function of real materials, including the ide- 
alized Drude formula We shall mention only those 
properties which are important in the derivation of our 
basic results; the proofs of theorems are given in the 
above textbooks. The vacuum case e{uj) = 1 is excluded 
from the discussion. 

Due to the causal relation between the displacement D 
and the electric field E, the dielectric function of every 
medium can be expressed as 



e{uj) = 1 + / dTe"^^G{T) 



(22) 



where the function G(t) is finite for all values of r, in- 
cluding zero. In particular, G(r) tends to as r — s- oo 
for dielectrics and it tends to Ana {a is the conductivity) 
as r ^ oo for conductors. The relation (|^^ has several 
important consequences. 

Let us first consider the frequency u; to be purely real. 
It follows from ([^ that e*{uj) = e{—uj). Denoting e{uj) — 
e'{Lo) + ie"(a;), where both the real e'{Lo) and imaginary 
t" {lo) parts are real numbers, we thus have 



(23) 



For any real material medium with absorption it holds 



£"(w) > 0(< 0) for Lj > 0(< 0); 



(24) 



the sign of e'iuj) is not subject to any physical restriction. 

If u is complex, ui = uj' + ioj" , the representation ([22l) 
tells us that e(w) is an analytic function of w in the upper 
half-plane w" > 0. Apart from a possible pole at cj = 
(for conductors), the analyticity extends also to the real 
oj axis. 

The function e{uj) does not take real values at any finite 
point in the upper half-plane, except on the imaginary 
axis. We can deduce from (|^^ that for any complex uj 



it holds e*(a;) = e(— w*). For purely imaginary lo = icu" 
indeed we find 



eiiu") = e*{iiu") 



lme{Lu) — for to — iuo'' 



(25) 

Moreover, on the imaginary axis, e(a;) decreases mono- 
tonically from eo > 1 (for dielectrics) or from oo (for con- 
ductors) at = iO to 1 at w = ioo. Hence, in particular, 
has no zeros in the upper half-plane. 
With regard to the symmetries for a real the 
expansion of the dielectric function around the origin lo = 
reads: for a conductor with conductivity cr > 



+ c + 0{u), 



(26) 



where the sign of the constant c is not restricted; for a 
dielectric medium of static dielectric constant eo 



e{uj) = Eq + icuj + 0{uj'^), c > 0, 



(27) 



where the positive sign of the constant c is fixed by the 
physical requirements (f24|) . 

In the limit \uj\ oo, a Taylor series expansion of 
G(t) around r = 0+ in Eq. ((22)) implies that, for both 
conductors and dielectrics, 



Re [e{uj) - 1] = O 



Ime(w) = O 



1 



(28a) 
(28b) 



IV. STATIC CHARGE CORRELATIONS 

In the general case of the plane contact between two 
media of dielectric functions ei(w) and £2(1^) pictured in 
Fig.l, the quantum formula for the Fourier transform of 
the surface charge correlation function ((5]) in the long- 
wavelength limit g ^ was derived in Ref. [H . Its static 
< = version can be reexpressed as 



/3^qu(0,-z) 



Im/(w), 



(29) 



where the form of the function /(w) depends on the con- 
sidered (retarded or non-retarded) regime. In the re- 
tarded case, /(w) = /qu (w) is given by two equivalent 
representations 



r 



1 



/qu(^) = f^5(^) I \ ( \ , ( \ r \ 
^ 47r^ Ki(wje2(w) -t- K2(a')ei(ti'j 



ei(a;)e2(t^) 



1 



Ki{uj)e2{uj) - K2{uj)ei{uj) 
An^"'"' q^[ei{io) + £2(0;)] - i0-^ei{u)t2[uj)/c' 

' ' ^ (30) 



€l{uo) £2(0;) / ' 
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Im CO 




FIG. 2: The contour in the complex frequency plane for t = 0. 



where g{uj) is defined in (j8b|) and the (complex) inverse 
lengths Ki and K2, one for each of the half-space regions, 
are given by 



'«i,2(w) = - ^ei,2(w), 



ReKi,2(w)>0. (31) 



In the non-retarded case, /(w) = fqu'^\uj) is obtained 
from pO|) by setting the speed of light c ^ oo. Since 
Ki = K2 = q in this limit, we get 

^ 47r^ V'^i(^) £2(w) ei(w) + e2(^^) 

The derivation procedures outlined below are applicable 
to both retarded and non-retarded regimes and so, when- 
ever possible, we use the simplified notation f{uj) to cover 

both functions /qu (w) and fqu\to). 

For a real frequency w, the symmetry relation e*{uj) = 
e(— w) implies that K*{io) — k,(—uj). Consequently, 

Re/(a.) =Re/(-^), Im/(^) - -Im/(-L^). (33) 



As w — > 0, 
Im/(0) = 0, 



Re/(0) 



47r2 Vei(O) £2(0) ei(0) + e2(0) 



(34a) 
.(34b) 



If w is complex and |w| oo, with the aid of the asymp- 
totic relations (|28ap and (j28b[) we find that 



lim /(w) = 0. 

|a;|— >oo 



(35) 



It is complicated to calculate the correlation function 
directly from the formula (|29p by expressing the imag- 
inary part of /((jj), then integrating over uj and finally 
taking the q ^ limit. We shall find the value of the 
integral of interest in another way, by using integration 
techniques in the complex plane and the analytic prop- 
erties of dielectric functions, summarized in the previous 
section. 



We can change slightly the path of integration, writing 

/ ^/(c) = -^nf{0) +V r —f{u), (36) 

where Ci is the path following the real axis, except it goes 
around the origin w = in a small semicircle in complex 
upper half-plane whose radius e tends to zero (Fig. 2). 
The first term on the rhs of ([55)1 is the contribution of 
the negatively oriented semicircle around the origin, V 
denotes the Cauchy principal value avoiding the origin. 



V / —f{uj) = lim / —f{uj) + / — 

(37) 

It is easy to see from the symmetry relations (|33p and 
the equality (|34ap that 



duj 



duj 



-Im/(cj) 



(38) 



We can close the path Ci by a semicircle at infinity C2 
(see Fig. 2) along which the integral J^^ duj f (oj) / cu is 
because of the asymptotic relation ([55]) . Denoting the 
closed contour as C, C — Ci U C2, and applying the 
operation Im to both sides of the relation (pS)) . we arrive 
at 



duj 



Im/(w) = 7r/(0) +Im 



duj 



(39) 



The integral over the contour C can be evaluated by using 
the residue theorem at poles {tOj} of the function /(w) 
in the lj upper half-plane bounded by C, 



Im 



Jc ^ , 



Res{f,ujj 



(40) 



provided that Iies{f,ujj)/ujj is real (which will be the 
case); Res denotes the residue. The static correlation 
function ([29]) is expressible as 



/3S'qu(0,(7) = ^/(0) + 27r 



Res(/,Ljj) 



3 



(41) 



where the value of /(O) is real, given by (j34b|l . The 
original algebraic task thus reduces to the problem of 
searching for all poles of the function f{uj) in the uj upper 
half-plane. 

Both the retarded and non-retarded ([5^ versions 
of the /-function contain g{uj) defined in (j8bp . Since g{uj) 
can be expanded in uj as follows [l6| 



g{u;) = 1 + 



^2' 



6- 



27r 



(42) 



it has in the upper half-plane an infinite sequence of sim- 
ple poles at the imaginary Matsubara frequencies 



i^j, Res{g,ujj) 



(j = l,2,...). (43) 
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The non-retarded /-function ([32|) has no further poles 
in the upper half-plane since the dielectric functions ei(w) 
and €2(0;) do not take there real values at any finite point. 
The retarded /-function ([50]) might have some further 
poles at points ui satisfying the equation 



icq? 



ei(tj) e2(w) 



(44) 



Interestingly, in the case of the jellium in vacuum (with 
real dielectric functions), this is just the dispersion re- 
lation for the surface plasmons (polaritons); for a recent 
review, see 17]. We are interested in the long-wavelength 
limit g ^ 0. If g = 0, the only solution of Eq. (^1)) is 
w = 0; this point is not inside the contour C due to 
the presence of the semicircle around w = 0. Now, let 
us study how the solutions of Eq. (|44|) "glue off" from 
w = when q is infinitesimal, but not identically equal 
to zero. We need the small- w expansion of l/e(w). It 
follows from the expansions (|26p and (|27p that for both 
conductors and dielectrics we can write 



1 

e~0) 



—---ibui + 0{u?), 5 > 0, (45) 



where the material constant h, b — 1/ {Attct) for conduc- 
tors and b = c/ep for dielectrics, is always positive, except 
for vacuum when 5 = 0. For small q, Eq. (j44p thus ex- 
hibits two solutions 



A. Retarded regime 



In the retarded case (l30l). we have 



F^K^.q) 



-y 

[ei(»0)-e2(ze,-)]^ 



(48) 



We recall from Sec. Ill that the values of the dielec- 
tric functions ei,2(jCj)i ^'^d consequently of the inverse 
lengths Ki,2(«Cj) (|3T|) . are real. We are interested in the 

1/2 

limit q ^ for which Ki.2(i^j) = Cj^i 2 Since cx j 

and, according to ()28a|) . e(jCj) — 1 = 0[l/p), the sum in 

(j48|) converges. This means that the function Fqn{Q,q), 
being of the order 0{q^), becomes negligible in compari- 
son with the first term in (|47ap when q ^ Q. In view of 
the representation ([6]), we find the static prefactor asso- 
ciated with the asymptotic decay to be 



1 



1 



ei(0) £2(0) ei(0) + e2(0) 



(49) 



Since this expression does not depend on the temper- 
ature and h, its classical f3h limit is the same, i.e. 



1 

8^ 



1 



1 



ei(0) £2(0) ei(0) + e2(0) 



(50) 



B. Non-retarded regime 



1 



ei(0) £2(0) 



1/2 



(46) 



Since 61 + 62 > 0, the two poles uj± move, as q increases 
from to a small positive number, from lu — to the 
lower u! half-plane, i.e. outside of the region enclosed by 
the C-contour. We conclude that, in both retarded and 
non-retarded regimes, only the poles on the imaginary 
axis at the Matsubara frequencies ((43|l contribute to the 
static correlation function ([4T|l . which thus becomes ex- 
pressible as follows 



/35qu(0,g) = 



FiO,q) 



An V£i(0) 
+F{0,q), 



£2(0) £i(0) + e2(0) 



(47a) 
(47b) 



The first term on the rhs of Eq. (|47aP is independent of 
(3h and c, the explicit form of the (static) function F{0, q) 
depends on the considered (retarded or non-retarded) 
regime. 



In the non-retarded case (|32p. we have 



27r^Vei(*0) e2(i^j) 



ei(«'?j) + £2(j^j) 



(51) 



It is evident from the asymptotic behavior £{i£,j) — 1 = 
0{l/j^) that the sum in (|5ip converges. The function 



Fqu''iO,q) is of the order 0{q) and its contribution to 
PS^\o, q) in (|47ap is nonzero in the limit q ^ 0. 

The explicit evaluation of the infinite sum over the 
Matsubara frequencies in ([5T|) is, in general, very compli- 
cated. As a check of the presented formalism, we recon- 
sider the previously studied case of the jellium in contact 
with vacuum, i.e. 



£l(cc>) = l-^, e2{Lu) = l. 



(52) 



Inserting these dielectric functions into ((5T|) and using 
the analog of the summation formula (I42p . 



E 



1 



{a coth a — 1) 



(53) 
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for a = f3fkjp/2 and phuJs/'2 {wg = ujp/\/2), we obtain 



Im CO 



) " ffH) - 1] 



(54) 



with g{uj) defined in (jSbjl . Adding to this function the 
first term g/(47r) in (|47ap . we reproduce the t = case 
of the previous result ((5a)) . 

In the classical (high-temperature) limit /3ft ^ 0, each 
of the frequencies tends to infinity, the corre- 

sponding terms in the summation over j in Eq. (|5ip 

vanish and so FqS'^'' (0, q) —i- 0. We are left with only the 
contribution identical to the retarded classical result (l50ll , 



i.e. /i^"'^^(0) = ft^[''(0), as it should be. For the configura- 
tion of a conductor, ei(0) = oo, in contact with vacuum, 
£2(0) = 1, we recover the universal result ([7]). For the 
configuration of a dielectric, ei(0) = eo, in contact with 
vacuum, €2(0) = 1, we recover our previous result (j2ip . 
The classical static prefactor /ici(O) to the l/B? asymp- 
totic decay is nonzero for an arbitrary configuration of 
two distinct media, except for the special case of two 
conductors. In that special case, the asymptotic decay is 
of a short-ranged type. 



V. TIME-DEPENDENT CHARGE 
CORRELATIONS 

For an arbitrary time difference t > between two 
points on the interface, the quantum formula for the 
Fourier transform of the surface charge correlation func- 
tion JSI) reads 



(55a) 



^e-^-*Im/(c.), 
— cos((jjt)Ini /(tj), (55b) 



where the retarded form of /{tu) is given in Eq. ([30| and 
the non- retarded one in Eq. ([32]) . 

Substituting f{uj) by e^'^'^* f {io) we can follow, in prin- 
ciple, the procedure outlined between Eqs. and ([55]) 
of the previous section to extend the integration over real 
frequencies to an integration over the contour C in Fig. 
2. The trouble is that on the Ihs of ([55)1 we end up with 
the integration over uj oi multiplied by 

Ime"*'^*/(cj) = cos(tjt)Im/(tj) +sin(cjt)Re/(w). (56) 



Comparing this expression with the representation (I55b|) 
we see that only the first term is needed. Within the 
C-contour formalism, we did not find a way how to get 
off the second (unwanted) term. 

Our new strategy is based on a transition from real 
times i > to imaginary times. In order to keep the 
function e~"^*/(tL')/w integrable over real lo, we make 
the substitutions t ^ —ir for u; > and t ^ ir for 
w < (r > 0), transforming in this way e~*"*/(^) ^ 
e~'"''^/(aj). In the w upper half-plane, the contour C in 




Re CO 



FIG. 3: Two contours in the complex frequency plane for 
t 0; u)j — i^j (j = 1,2,...) are the Matsubara frequencies 



Fig. 2 will be replaced by two contours C+ and C_ in the 
quarter-spaces Rew > and Reoj < 0, respectively (see 
Fig. 3). The contours are constructed in such a way that 
one avoids the singularities (simple poles) of f{uj)/uj, the 
one at w = and the infinite sequence of poles at the 
imaginary Matsubara frequencies The contour 

C-f is directed along the real axis from to cx), except for 
an infinitesimal quarter-circle around the origin lu = 0, 
continues by a quarter-circle at infinity and returns to 
the origin neighborhood along the imaginary axis, avoid- 
ing by infinitesimal semicircles the Matsubara frequen- 
cies The contour C_ is directed along the real axis 
from —00 to 0, except for an infinitesimal quarter-circle 
around the origin w = 0, continues along the imaginary 
axis, avoiding by infinitesimal semicircles the Matsubara 
frequencies {i^j}, and returns to the starting point by a 
quarter-circle at infinity. In close analogy with the proce- 
dure outlined between Eqs. (|36p and ([39)) . we can derive 
the integral equality 



iTTf{0)+V 



C- 



OJ 



cLu 



(57) 



Here, we have used two facts: the contributions com- 
ing from the straight-line fragments of the paths Cj^. and 
C- between two neighboring frequencies (i^j, i^j+i) can- 
cel exactly with one another due to the opposite direc- 
tion of the integration and the contribution of the C+- 
semicircle (C_-semicircle) around the pole i^j is equal to 
— z7rRes(/, multiplied by the corresponding time- 

dependent factor e-*«^^ (e*«^^). The Ihs of Eq. ^ is 
equal to zero since there are no poles of the integrated 
functions inside the contours C_|_ and C_. Since the ra- 
tio Res(/, is a real number, taking the imaginary 
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part of Eq. ((57|) leads to 



we obtain 



/OO 7 
— e-H-Im/M=^/(0)+F(r,g), (58a) 
-OO ^ 



where 



(58b) 



The function F{t, q) represents the (imaginary) time- 
generahzation of the static F{0, q) (j47bp . In the retarded 
case, the generahzation of the static formula reads 



F^lli.r,q) 



''''''' 27r ^ Ki(iCi)e2(i^j) + K2(i0)ei(«0) 

^ h(.0) -..(<,)]- 
ei(«0)e2(«?j) 

In the non-retarded case, the generalization of the static 
formula (ISD) reads 



1 



1 



ei(iCi) + £2(^6) 



cos(^jr). (60) 



We have shown in Sec. IV that the series determining 
the functions in (|59p and ([50)1 are convergent for r = 0. 
The presence of the oscillating factor cos(^jr) in the se- 
ries for r > even improves their convergence property. 
Let us "pretend" that we have found the explicit form of 
the function F{t, q) (|58b[) . the type of the regime is irrele- 
vant. To express the time-dependent correlation function 
(j55b[) . we have to return from imaginary to real times by 
considering the substitution r — s- it in Eq. (|58ap . Since 
the integral on the Ihs of (|58ap is finite for all complex 
T with Rer > 0, there must exist a well-behaved ana- 
lytic continuation F{it,q) of the function F{T,q). Con- 
sequently, 



PSquit,q) = — 



1 



1 



47r Vei(O) £2(0) 
+ReF{it,q). 



ei(0)+e2(0) 



(61) 



For t = 0, the operation Re on F{0, q) become superflu- 
ous and we recover the static result (|47al) . 

To check that the formalism works correctly, we recon- 
sider the non-retarded regime for the contact between the 
jellium and vacuum. Inserting the corresponding dielec- 
tric constants ([5^ into ([SO)) and using the summation 
formula UM 



{a/Try 



(j-UJT 

\ a 



cosh(a — lot) 
a — 1 



sinha 



(62) 



cosh[{(3Hujs/2) — cUst] 



- 1 . (63) 



si-nh{PhuJs/2) 
Phup cosh[{/3hujp/2) — oJpT] 
~2 smh{phujp/2) 



The analytic continuation of this function from r to it 
is well defined. Using the relation Re cosh(a — iojt) = 
coshacos(a;t) valid for real a and ojt, we get 

Rei^^-H*^,?) - ^[2giu;s)cos{u;,t) 

-g{ujp) cos(wpi) - 1] (64) 

with g{uj) defined in (j8bp . Adding to this result the first 
term q/{A'K) in (|6ip . we reproduce the previous time- 
dependent result ([5a|) . 

In the retarded case, the function Fqu (t, q) in ([55)1 is of 
the order 0{q'^). The same property holds for its analytic 
continuation Fqu{it,q). Eq. ([61]) then tells us that, in 
the limit g — > 0, 



/3^W(t,, 



4^ Vei(O) £2(0) £i(0) + £2(0) 



(65) 

In view of the representation ([6]), the quantum time- 
dependent prefactor to the asymptotic decay takes, for 
any temperature, its static classical form 



1 



/iWm = L [ ^_ + 

'^"^^ 8^2 |^£l(0) ^ £2(0) £l(0)+£2(0) 



(66) 



This result holds for all possible media combinations if 
one takes £(0) ioo for conductors, £(0) = £o > 1 for 
dielectrics and £(0) = 1 for vacuum. 



VI. CONCLUSION 

Although the present work is rather technical, its main 
result ([^ (or, equivalently, ((Ml)) is of physical inter- 
est. It represents the generalization of the analogous re- 
sult, obtained in Refs. [ll, Q and valid exclusively for 
the special jellium model of conductors in contact with 
vacuum, to all possible media (conductor, dielectric, vac- 
uum) configurations. The derivation of the general result 
()65p . based on the integration in the complex plane and 
on the known analytic properties of dielectric functions 
in the frequency upper half-plane, is much simpler and 
more transparent than the one for the special case of the 
jellium in contact with vacuum [HQ. 

There is still an open question: which are physical 
reasons that the inclusion of retardation effects causes 
the asymptotic decay of time-dependent quantum surface 
charge correlations to take its static classical form, inde- 
pendent of h and c? For the static quantum case t = 0, 
a possible explanation might be that the surface charge 
correlation function depends on the only dimensionless 
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quantity phc/R constructed from universal constant and 
the distance. If this is the case, the large-distance asymp- 
totics i? cx) is equivalent to the classical limit (3h ^ 0. 
A verification whether this claim is true or not is left for 
the future. The case t ^ will be still an open problem. 
It might be also possible to verify our results by using a 
fully microscopic approach. 

We believe that Rytov's fluctuational electrodynamics 
and the techniques presented are applicable also to other 
important phenomena associated with the presence of a 
surface between distinct media. One of such problems 
can be the evaluation of the shape-dependent dielectric 
susceptibility tensor [§, [H, [l^ for quantum Coulomb flu- 
ids, with and without retardation effects. Another inter- 



esting topic is the large-distance behavior of the current- 
current correlation function near an interface between 
two media. 
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